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Abstract. We present simulations of Far Infrared (FIR) 
emission by dust in spiral galaxies, based on the Monte 



Carlo radiative transfer code of Bianchi, Ferrara & Giova- 
nardi (1996| ). The radiative transfer is carried out at sev- 



eral wavelength in the Ultraviolet (UV), optical and Near 
Infrared (NIR), to cover the range of the stellar Spectral 
Energy Distribution (SED). Together with the images of 



distribution, the extinction of starlight may affect to dif- 
ferent degrees our knowledge of the objects harboring dust 
and of the distant universe in the background. However, 
deriving the amount of dust from the extinction features 
in spiral galaxies is not easy, because of the uncertainty in 
the relative distributions of dust and stars. Indeed, realis- 



the galactic model, a map of the energy absorbed by dust 



is produced. Using Galactic dust properties, the spatial 
distribution of dust temperature is derived under the as- 
sumption of thermal equilibrium. A correction is applied 
for non-equilibrium emission in the Mid Infrared (MIR). 
Images of dust emission can then be produced at any wave- 
length in the FIR. 

We show the application of the model to the spiral 
galaxy NGC 6946. The observed stellar SED is used as 
input and models are produced for different star-dust ge- 
ometries. It is found that only optically thick dust disks 
can reproduce the observed amount of FIR radiation. 
However, it is not possible to reproduce the large FIR 
scalelength suggested by recent observation of spirals at 
200^m, even when the scalelength of the dust disk is larger 
than that for stars. Optically thin models have ratios of 
optical/FIR scalelengths closer to the 200//m observations, 
but with smaller absolute scalelengths than optically thick 
cases. The modelled temperature distributions are com- 



tic models are needed to avoid equivocal results (Disney 
ct al. 1989| ) . Recently, Xilouris ct al. (1999) have analised 



a sample of seven edge-on spirals by fitting optical and 
NIR images with a complex radiative transfer model in- 
clusive of scattering. They derived a small opacity, with 
their exponential dust disks having a mean face-on optical 
depth in the B-band tb = 0.8. Unfortunately, the method 
can only be applied successfully to a limited numer of 
cases, i.e. edge-on galaxies, where it is possible to analyze 
the vertical surface-brightness distribution and the dust 
extinction is maximised by the inclination. 

Further constrains on the structure and opacity of 
dust disks can be obtained by studying dust emission. As- 
suming that all the stellar radiation absorbed by dust is 
reemitted in the MIR and (mainly) in the FIR, the amount 
of extinction in a galaxy can be simply derived by com- 
paring the total dust emission with the stellar one (the 
energy balance technique) . The extinction is then used, to- 
gether with a radiative transfer model, to study the dust 



patiblc with observations of the Galaxy and other spirals. 1998a). However, the dust emission is treated as a bulk 



disk properties (Evans 1992; Xu & Buat 1995; Trewhella 



We finally discuss the approximations of the model and 
the impact of a clumpy stellar and dust structure on the 
FIR simulations. 

Key words: Radiative transfer - dust, extinction - 
Galaxies: individual: NGC 6946 - Galaxies: spiral - Galax- 
ies: structure - Infrared: galaxies 



1. Introduction 

Dark lanes in edge-on spirals clearly reveal the presence of 
dust in galactic disks. Depending on the dust amount and 



Send offprint requests to: bianchi@mpia-hd.mpg.de 



and generally only the disk optical depth is derived. The 
parameters for the geometrical distributions of dust and 
stars have to be assumed a priori. 

Three-dimensional models of FIR emission from spiral 
galaxies can help to retrieve information about the rel- 
ative star-dust geometry from the SED of dust emission 
and the spatial distribution of the radiation, when avail- 
able. However, an additional factor is introduced in such 
models. Not only it is necessary to assume a dust distri- 
bution, but also a knowledge of the variation of the ISRF 
heating throughout the disk is needed. The ISRF could be 
in principle derived from the observed profiles of stellar 
emission, but a complicate de-projection is needed involv- 
ing assumption on the dust distribution. This approach is 
adopted by Walterbos fc Greenawalt (1996 ). They com- 
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pute fluxes and surface brightness gradients at 60/im and 
100/im for a sample of 20 spirals. The diffuse ISRF is de- 
rived from B band profiles and the dust distribution from 
the HI column density. The problem of defining an ISRF 



not yet scientifically validated (Klaas & Richards 2000) 



We will use the model presented here to test if indeed 
the observations are compatible with large disks. Other 



works seems to support the same hypothesis. Davies et al 



can be avoided by assuming a temperature gradient for (1997) could model the Galactic FIR surface brightness 



dust grains. Davies et al. (1997) modelled the 140/im and 
240/im Galactic emission observed by DIRBE, by adopt- 
ing a gradient that matches the longitudinal variation of 
temperature observed in the data themselves. 

Otherwise, the dust heating by the ISRF can be de- 
rived in a self-consistent way from a radiative transfer 
model of the stellar radiation in the dusty environment. 



Sauty et al. (1998) produced a Monte Carlo radiative 



transfer model for the clumpy ISM of the spiral galaxy 
NGC 6946, although a self-consistent ISRF is derived only 
for the UV radiation, relying on optical profiles for ra- 
diation at longer wavelelngth. Dust heating in a clumpy 
galactic ISM is also included in the photometric evolution 



model of Silva et al. (1998). They carry out the radiative 
transfer in a proper way for stellar emission embedded 
in the clumps and use an approximation for the smooth 
medium. In the emission model of |Wolf et al. (1998| ) the 
radiative transfer is computed with the Monte Carlo tech- 
nique for the whole spectral range, both in the smooth 
medium and clumps. However, the adopted geometry is 
typical of star-formation environment. 

In this paper, we present a self-consistent galactic 
model for the FIR emission. The model is based on the 
Bianchi, Ferrara fc Giovanardi (1996j , hereafter, BFG) 



Monte Carlo code for the radiative transfer (complete 
with scattering) in dusty spiral galaxies. A map of the 
energy absorbed by dust is derived from the output of 
the radiative transfer for stellar radiation at several wave- 
length. The grain temperature along the dust distribution 
is therefore computed from the radiative transfer itself. 
Luminosity, spectrum and surface brightness distribution 
of the FIR radiation can be retrieved, for any wavelength. 
The main advantage with respect to other models of FIR 
emission lies in the relative simplicity. We have limited 
the number of parameters involved to a minimum, adopt- 
ing, whenever possible, empyrical data. This will allow 
an isolation of the effect of the dust distribution on the 
FIR emission. We will try to explain the characteristic of 
FIR emission with diffuse dust heated by a diffuse ISRF. 



Clump ing of dust and embedded stellar emission ( Bianchi 
et al. 2 000b| ) are not yet included in this work. 



Recent FIR observations suggest large scalelengths for 
dust emission. In a sample of seven spirals resolved by the 



ISOP HOT instrument aboard the ISO satellite (Kcsslcr 
et al. 1 |996); Lcmkc et al. 1996[ ), |Alton et al. (1998cj ) found 



that the 200 /im scalelengths are 70% larger than those 
in the B-band. The large ratio may be a result of a dust 
disk larger than the stellar distribution. However, the re- 
sult may also be affected by the transient effects of the 
ISOPHOT detectors. In fact, th e P32 mapping mode with 
which the Alton et al. (1998c ) images were obtained is 



only with a dust disk with radial and vertical scalelengths 
larger than the stellar by a factor 1.5 and 2.0, respectively. 
Nelson, Zaritsky fc Cutri (1998| ) combined 100 fim IRAS 
images of galaxies with similar angular size and concluded 
that dust is present w ell beyond the observed stellar disk. 
Xilouris et al. (1999| ) derived intrinsic dust scalelengths 
larger than the stellar one, by a mean factor 1.4 in the 
seven edge-on spiral sample. If confirmed, such disks may 
have a large impact on the observations of the distant 
universe, because of their larger cross section to radiation 
from backgroud objects. 

The paper is structured as follows: Sect. || describes 
the radiative transfer code, with a particular attention 
to the chosen dust and stellar distributions and to the 
modification with respect to the orig inal BFG code. Sect.| 
and ^ show how the map of absorbed energy derived from 
the radiative transfer model can be converted into a map 
of the dust temperature distribution and hence, of FIR 
emission. We then show an application of the model to 
the spiral galaxy NGC 6946 with the intent of analysing 
the behaviour of our model and test the influence of the 
various parameters on the FIR emission. The observed 
properties of NGC 6946 used in the model are described in 
Sect. ||. The results are presented in Sect. || and discussed 
in Sect. [?]. A summary is given in the last section. 



2. The radiative transfer model 

The original BFG radiative transfer code carries out the 
radiative transfer for typical galactic geometries. Simula- 
tion of optical surface brightness, as well as light polari- 
sation, are produced. Dust properties are computed from 
the Draine fc Lee (1984 ) dust model, using Mie's theory 
for spherical grains. In this paper we use a simplified ver- 
sion of the code: the treatment of polarization has been 
omitted; empirical dust properties and scattering phase 
functions are used. The same approach has recently been 
adopted in Ferrara et al. (1999) and Bianchi et al. (2000b). 

In this Section, we describe the stellar and dust distri- 
butions adopted in the galaxy model, and the chosen dust 
extinction and scattering properties. A brief description is 
then given of the radiative transfer code. 



2.1. The stellar disk 

As in BFG, the new version of the code allows the use of 
both disks and spheroidal bulges to describe the spatial 



distribution of stars (Bianchi 1999). However, to limit the 



number of free parameters in the modelling, we adopt here 
a single stellar disk. The chosen disk is exponential along 
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Martinet 1994), we derive a*/(/3*) = 14.4. Similar values 



density distribution is thus described by 


for a*/(/3*) can be found in the literature ( 


van der Kruit 


P - Po- 




& Searle 1982 


: de Grijs & van der Kruit 199£; Xilouris 


!Ap( r/a*- | i | //?*)■ 


W etal. 199£). 







where r and z are the galactocentric distance and the 
height above the galactic plane, respectively, and a+ and 
/3* the relative scalelengths. While there is a consensus for 



The stellar disk is truncated at 6a*. Fits of star counts 
f or faint Galactic sou r ces suggest a similar trunc a tion 
(IWainscoat et al. 1992| ; |Robin, Creze fc Mohan 1992j ; [RaT 



the ra dial exponential behaviour (|dc Vaucouleurs 1959|-| phy et al. 1996|) . The truncation along the vertical direc 



Freeme n 1970 ), a number of expressions fo r the vertical tion is at 6/3+ 



distrib ution ha ve been used in the literature, van der Kruit 
& Sear |.c (1981 ) use a sech 2 , the solution for a self gravitat 



ing isothermal sheet, while van der Kruit (1988 ) propose 
a sech, consistent with measures of stellar velocity disper- 
sion out of the Galactic plane. The expo nential adopted 
here ( Wainscoat, Freeman fc Hyland 1989| ) has the advan- 
tage of mathematical simplicity, but has no firm physical 
justification. Analysing the vertical profiles of a sample of 
24 edg e-on galaxies in the relatively du st-free K-band, [ie] 
Grijs, 3 eletier fc van der Kruit (1997 ) find that a distri- 



bution with a peak intermediate between the exp and the 



sech fits better the central peak (see also de Grijs & van 



der Kruit 1996). Since a small inclination from the pure 
edge-on case can produce a less sharp profile, these results 
are consistent with an exponential distribution. 



The radiative transfer model of the galaxy (Sect. 2.4) is 
run for several wavelength bands. The radial scalclcngth of 
the stellar distribution, a*, is assumed to be the same for 
all of these bands. Therefore, any observed colour gradient 
along the disk of a spiral is interpreted as due to differ- 
ential extinction with A, rather than being a reflection of 
different stellar components. The wavelength dependence 
of the intrinsic stellar radial scalelength is uncertain, be- 
cause a knowle dge of the dust distr ibution is needed for its 
determination. Peletier et al. (1995 ) analyzed the variation 
of the ratio of B and K-band scalelengths with inclination 
in a sample of 37 Sb-Sc galaxies and concluded that the 
observed gradient is due to extinction. On the other hand, 
do Jong (1996 ) compared colour-colour plots of a sample of 
86 face-on galaxies with a Monte Carlo radiative transfer 
model and finds that reasonable dust models cannot be re- 
sponsible for the observed colour gradients. However, such 
a conclusion does not seem to include models with dust 
disks more extended that the stellar, that is one of the 



cases studied here. The fits of edge-on galaxies of Xilouris 



et al. (1999) show a slight variation of the intrinsic scale- 
length with wavelength. 

Observations on the Galaxy show that different stellar 
populations have different vertical scalelengths. For the 
sake of simplicity, a single mean value for the vertical 
scalelength has been adopted here, for each wavelength. 



Wainscoat et al. (1992) provide a table of vertical scale- 
lengths, B, V, J, H, K-band absolute magnitudes and lo- 
cal number density for the main Galactic stellar sources 
A mean value has been computed averaging over the to- 
tal disk luminosity inte grated from the B to the K band 
Assuming a+ = 3 kpc ( Kent, Dame fe Fazio 1991; Fux fc 



2.2. The dust disk 

The parameters for the dust disk are far more uncertain 
than those for the stellar disk. Usually the same functional 
form is used as for the stellar distribution, with indepen- 



dent 


scalelengths (Kylafis & Bahcall 1987 


; Byun et al. 


1994 


; Bianchi et al. 1996; Xilouris et al. 1999 


)■ 



A very good correlation is usually found between dust 
and gas in spirals, especially with the molecular compo- 



ncnt which is dominant in the inner galaxy (Devereux & 


Young 1990|; Xu et al. 1997 


; Bianchi et al. 1998 


; Alton 


3t al. 1998a; Bianchi et al. 2000a 


). It is therefore reason- 



able to use the gas distribution as a tracer of dust. 

In luminous, face-on, late-type spirals, H2 peaks in the 
centre and falls off monotonically with increasing galacto- 
centric distance, while the HI distribution shows a central 
depression and a nearly c onstant surface density across 
the rest of the optical disk ( [Young fe Scovillc 1991 ) . Some 



early type galaxies show the same behaviour, although a 
good fraction of them presents a central depression and a 
flatter ring distribution for the molecular gas. Indeed, flat 
or ring-like dust distributions have been derived from FIR 



and sub-mm observations ( 


Xu & Helou 19961; Sodrosk 


|et al. 1997 


; Bianchi et al. 1998 


)■ 



Nevertheless, the gas distribution of most galaxies is 
dominated by a centrally peaked molecular gas compo- 
nent. This is the case of the late- type spiral NGC 6946 



(Tacconi & Young 1986), whose dust distribution is the 
main concern of this paper (Sect. |^). Therefore, in this 
paper dust is assumed to be distributed in a smooth ra- 
dial and verti cal exponential disk, similar to the stellar 
one (Sect. 2T). The number density of dust grains can be 
written as 



n(r,z) = n exp(-r/a rf - | z \ //3 d ), 



(2) 



with n,Q the central number density. The dust scalelengths 
aj, and fid can be selected independently from the analo- 
gous stellar parameters. It is usually assumed in radiative 



transfer models that ad ~ a* and (3d ~ 0.5/3* 1 


Byun et al. 


1994 


; Bianchi et al. 1996 


; Davies et al. 1997 


and refer- 



ences therein) . We will refer to this as the standard model. 
The choice of the vertical scalelengths is mainly dictated 
by the impossibility of simulating the extinction lanes in 
edge-o n galaxies with a du st distribution thicker than the 
stars ( Xilouris et al. 1999 ). The assumption of a similar 
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radial ucalclcngth for stars and dust has no firm justifica- 
tion. As already said in the introduction, recent analysis 
of extinction of starlight and FIR emission suggests more 
extended dust distributions. 

The central number density of dust no is normalised 
from the chosen optical depth of the dust disk. For a cen- 
tral face-on optical depth in the V-band, -ry, no is given 

by 



n 



TV 



'2P d a ext (V) , 



(3) 



where o~ ex t(V) is the extinction cross section. The optical 
depth at any wavelength can then be computed by inte- 
grating the absorption coefficient 



k\(r,z) = n(r, z)a ext (X) 



(4) 



along any path through the dust distribution. The integral 
involves knowledge of the ratio croxt(A)/cr cx t(V^), which is 
found from the assumed extinction law. As for the stellar 
disk, the dust disk has been truncated at 6 scalelengths 
both in the vertical and radial direction. 

In section ^various values of ad/ a*, (3d/ (3^ and ry will 
be explored. 

2.3. Dust extinction 

The extinction law adopted here is that typical of dif- 
fuse Galactic dust, with a strong bump at 2175A. Re- 
gions of high star formation and starburst galaxies usually 



have a weaker bump and a steeper far-UV rise ( 


Whittet 


1992; Calzetti, Kinney & Storchi-Bcrgmann 1994 


; Gor- 


don, Calzetti & Witt 1997). Extinction laws in the optical 



show smaller differences. From radiative transfer models 
of seven edge-on galaxies, Xilouris et al. (1999] ) derive an 



extinction law similar to the Galactic one longward of the 
U-band. 

The extinction law used is tabulated in Tab. ([!]), for 
17 bands in the spectral range of stellar emission. The 
bands from UV1 to K have essentially the same spectral 



coverage as the homonymous described by Gordon et al 



(1997 ), from which the extinction law is taken. Two bands 
(namely EUV and LMN) have been added to extend the 
spectral coverage in the ultraviolet up to the ionization 
limit and in the near infrared. The extinction law for these 



two bands has bee n taken from Whittet (1992 ) and Rieke 
fc Lebofsky (1985 ), re spectively. 



We have used the Henyey fc Greenstein (1941 ) phase 
function for the scattering. This is given by 



1 



1 - 9 



2(l+g 2 - 2gcos9fl 2 '' 



(5) 



with 9 the angle of scattering. The asymmetry parameter 
<7, as well as the dust albedo u>, are derived from models 
of dust scattering by reflection nebulae in the Milky Way. 
The adopted g and u) are shown in Tab. (f[ ). Again, the 



values for the bands from UV1 to K are from Gordon et al 



(1997| ). For the EUV band we have used |Witt et al. (1993Q 
data at 1000A while the values for the LMN band have 
been assumed equal to those in the K band. 

2.4- The Monte Carlo code 

We give here a brief description of the radiative transfer 
code, referring the interested reader to the BFG paper for 
more details. 

Adhering to the Monte Carlo method, the code follows 
the life of each energy unit (a photon) through scattering 
and absorption, until the radiation is able to escape the 
dusty medium. The code is monochromatic, since the op- 
tical properties of dust must be specified for a particular 
wavelength. In principle, the geometrical distributions of 
stars may depend on A too. For a given star-dust geometry 
and central face-on optical depth ry, a radiative transfer 
simulation is produced for each of the 17 bands of Sect. |2.3| . 
Typically, a simulation consists of 10 7 photons. 

The main steps of the computation scheme are: 

Emission: the position of a photon in the 3-D space is 
derived according to the stellar distributions described 
in Sect. 2.1. Photons are emitted isotropically and with 
unit energy. 

Calculation of optical depth: The absorption coefficient 
k\(r, z) (Eqn. [|) is integrated from the emission posi- 
tion along the photon travelling direction to derive the 
total optical depth tt through the dust distribution. 
A fraction e~ TT of all the energy travelling in that di- 
rection propagates through the dust. With the Monte 
Carlo method it is then possible to extract the optical 
depth r at which the photon impinges on a dust grain. 
This optical depth can be computed inverting 



d£, = R, 



(6) 



where R is a random number in the range [0,1]. If the 
derived r is smaller than tt, the photon suffer extinc- 
tion, otherwise it escapes the dusty medium. This pro- 
cess is quite inefficient when the optical depth of the 
dust distribution is small, most of the photons leav- 
ing the dust distribution unaffected. To overcome this 



problem, the forced scattering method is used (Cash 



well fc Everett 195£; Witt 1977): essentially, a fraction 



e~ TT of the photon energy is unextinguished and the 
remaining 1 — e~ TT is forced to scatter. When the op- 
tical depth is small (tt < 10~ 4 ) or the photon path 
is free of dust, the photon escapes the cycle. Once r 
is known, the integral of k\(r,z) is inverted to derive 
the geometrical position of the interaction between the 
photon and the dust grain. 
Scattering and Absorption: a fraction of the photon en- 
ergy, given by the albedo u>, is scattered, while the 
remaining (1 — u>) is absorbed. The scattering polar 
angle 9, i.e. the angle between the original photon 
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band 


A; < A [A]< A f 


0- C xt(A)/cr ox t(V) 




9 


-Fmir(A) 


EUV 


912 


1125 


4.34 


0.42 


0.75 


0.72 


UV1 


1125 


1375 


3.11 


0.60 


0.75 


0.54 


UV2 


1375 


1655 


2.63 


0.67 


0.75 


0.44 


UV3 


1655 


1900 


2.50 


0.65 


0.73 


0.47 


UV4 


1900 


2090 


2.78 


0.55 


0.72 


0.58 


UV5 


2090 


2340 


3.12 


0.46 


0.71 


0.65 


UV6 


2340 


2620 


2.35 


0.56 


0.70 


0.54 


UV7 


2620 


3230 


2.00 


0.61 


0.69 


0.43 


U 


3230 


3930 


1.52 


0.63 


0.65 


0.36 


B 


3930 


4979 


1.32 


0.61 


0.63 


0.28 


V 


4979 


5878 


1.00 


0.59 


0.61 


0.23 


R 


5878 


7350 


0.76 


0.57 


0.57 


0.18 


I 


7350 


9500 


0.48 


0.55 


0.53 


0.14 


J 


9500 


14000 


0.28 


0.53 


0.47 


0.12 


H 


14000 


19000 


0.167 


0.51 


0.45 


0.12 


K 


19000 


25000 


0.095 


0.50 


0.43 


0.12 


LMN 


25000 


37500 


0.04 


0.50 


0.43 


0.12 



Table 1. Milky Way extinction law (c e xt(A)/(7ext(10)> albedo lo and asymmetry parameter g and wavelength range 
of the 17 bands used in the models. The last column gives the value of the MIR correction -Fmir (Eqn. ||). 



path and the scattered direction, is computed using the 
|Hcnycy fc Grccnstcin (1941 ) scattering phase function 
(Eqn. ^|), inverting 



9') sin9' = i?, 



(7) 



with R a random number. The inversion of Eqn. (@) is 
given by the analytical formula 



1 

~9 



1 



(1 



Jl\2 



(l+g(l-2R)f 



(8) 



No preferential direction perpendicular to the original 
photon path is assumed, as for scattering by spherical 
grains. The azimuthal angle is thus extracted randomly 
in the range [0, 2ir}. 

The amount of energy absorbed by the dust grain is 
stored as a function of the galactocentric distance and 
the height above the plane, using the two model sym- 
metries, i.e. the symmetry around the vertical axis and 
the symmetry about the galactic plane, to improve the 
signal-to-noise. Maps of absorbed energy were not pro- 
duced by the original BFG code. 
Exit conditions: the last two steps are then repeated, us- 
ing the new direction of the scattered photon, the coor- 
dinates of the scattering point and the energy reduced 
by absorption, until the energy of the photon falls be- 
low a threshold value (10~ 4 of the initial energy) or 
until the exit conditions on r are verified. 

After the exit conditions are satisfied, the photon is char- 
acterised by the last scattering point, its travelling direc- 
tion and its energy. To reduce the computational time, the 
two model symmetries are exploited to produce a total of 
4 photons from each one. Photons are then classified ac- 
cording to the angle between the last direction and the 



symmetry axis. An image is produced collecting all the 
photons that fall in a solid angle band of width 47r/15 
(BFG) and mean polar angle equal to the chosen model 
inclination. 

Finally, all the photons in an angle band are projected 
into the plane of the sky according to their point of last 
scattering. Maps of 201x201 pixels are produced, cover- 
ing an area of 12x12 stellar radial scalclcngths around the 
centre of the galaxy. Maps of absorbed energy cover 6 dust 
scalelengths in the radial direction and in the positive ver- 
tical direction in 101x101 pixels. 

2. 5. Normalisation of the radiative transfer output 

Each set of 17 radiative transfer simulations is then nor- 
malised according to a chosen SED for the stellar radia- 
tion. Two kinds of normalization are available: in the first, 
the input normalization, a SED is chosen for the intrinsic, 
unextinguished, stellar energy emission. The SED is inte- 
grated over the bands limits and the total energy emitted 
in a band is used to scale each monochromatic simulation. 
This is suitable, for instance, to predict the FIR emission 
associated with a synthetic galactic SED. 

The second normalisation mode, the output normalisa- 
tion, is instead based on the observed SED of a galaxy. An 
observed SED is constructed from observations in several 
wavelength and the model is scaled to produce images that 
have the same SED. The intrinsic, unextinguished, SED 
is then inferred from the radiative transfer model. This 
second mode, suitable for fitting the FIR emission, is the 
one used in this paper (Sect. ||). 

Together with the optical images, the absorbed energy 
maps are also normalised. Thus, the output of a model 
consists of a set of 17 images in the wavelength range 
of stellar emission and a set of 17 maps of the energy 
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absorbed by dust from stellar radiation emitted in each 
band. The images are produced in units of surface bright- 
ness, while the maps of absorbed energy are in units of 
energy per unit time per unit volume. The absorbed en- 
ergy maps can be added together to produce a map of the 
total energy absorbed by dust. However, for the purpose 
of modelling the FIR emission, a further step is needed 
before coadding. This is explained in the next Section. 



3. The MIR correction 

The energy absorbed from photons heats up dust grains 
and is re-emitted in the infrared, preferentially at A > 
10/im. We do not consider here the heating of dust grains 
by collisions with the interstellar gas: this process is nor- 
mally negligible and the dust temperatu re is almost en- 
tirely determined by radiative processes ( Bpitzer 1978 ). 

When the energy of the absorbed photon is small com- 
pared . to the internal energy of the grain, radiation io 



l/X [urn -1 ] 

Fig. 1. The contribution of very small grains and PAHsto 
the absorpt ion, -Fmir(A) (Eqn||), derived from the Desert 



emitte d at the thermal equilibrium. Thio io the caoo for 

larger grains, responsible for most of the FIR emission. For 
smaller grains, the absorption of a single high-energy pho- 
ton can substantially alter the internal energy. The grain 
undergoes temperature fluctuations of several degrees and 
cools by emission of radiation, mainly in the MIR range 
( Whittet 1992 ). The models of this paper are limited to 
thermal FIR emission only. Therefore, it is necessary to 
exclude from the total absorbed energy the fraction that 
goes into non-equilibrium heating. 

Desert, Boulanger & Puget (1990) derived an empir- 
ical dust model, from an analysis of the features in the 
Galactic extinction law and in the infrared emission. Three 
dust components were needed in the model: i) big grains 
(0.015/im < a < 0.11/im), responsible for the FIR emis- 
sion; ii) very small grains (0.0012/im < a < 0.015/im) and 
iii) PAHs, responsible for the MIR emission at A < 80/^m 
From the parameters and functional forms of Desert et al^ 
|(1990|) , we have derived the mean absorption cross-section 
of the model, er a bs(A), and for each of its three dust com- 
ponents^]. Of the light impinging on the grain mixture, 
a fraction of energy proportional to <7 a b s (A) is absorbed; 
therefore, the contribution of very small grains and PAHs 
to the absorption is given by 



et al. (1990| ) model 



After multiplying by (1 — -Fmir(A)), the 17 absorbed 
energy maps store the energy that is absorbed by large 
grains only. The final coadded map, W a bs (r, z) , contains 
the energy absorbed by dust per unit time and unit vol- 
ume, that goes into thermal equilibrium emission only. 
The derivation of the dust temperature and FIR emission 
from VFabs (r, z) is shown in the next Section. 



4. The dust emission model 

The power absorbed by a single dust grain can be derived 
by dividing W a bs(»", z) by the grain number density n(r, z). 
Because of the conservation of energy, this is equal to the 
power each grain radiates. Assuming that all dust grains 
in our model are spherical and with the same mean radius 
a, the thermal equilibrium imposes that 



Wa bs (r, z) 
n(r, z) 



Q em (\)irBx (T d (r,z))d\, 



(10) 



■Pmir(A) 



VSG 
abs 



(A) 



PAHs 
abs 



(A) 



Cabs (A) 



(9) 



Values of Fmir(A), the MIR correction, are given in 
Tab. (0) for each of the model bands and plotted in Fig. [j] 
as a function of 1/A. i 7 MiR(A) has a local maximum in the 
position of the 2175A bump of the e xtinction curve: ab- 
sorption by very small grains in the Desert et al. (1990 ) 
model is responsible for this feature. The rise in the Far- 
UV is due to PAHs. 



where T^(r,z) is the dust temperature distribution and 
Qem(A) the emission efficiency. By Kirchhof's law, Q om = 
Qabs = Cabs/™ 2 QWhittct 1992| ). 

For the exponential dust distribution of Eqn. (||) and 
the central number density of Eqn. (||), Eqn. ([!(]) can be 
rewritten as 



f3dW a bs(r, z) 



2T V exp(-r/a d - \ z \ / f3 d ) 



Qom(A) 
Qoxt(V0 



irB x (T d (r, z))d\, (11) 



Note that a cxt = o"abs + o" sca . In the Desert et al. (1990 ) 
model, very small grains and PAHs are pure absorbers (i.e. 

0"sca = 0). 



with Qext(^) = &ext(V) /ira 2 . The map of dust tempera- 
ture T<j(r, z) is derived by inverting Eqn. (pd|). 
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For the ratio Q e m (A) /Qext ( V) ; we use a value derived 
for Galactic dust ([Bianchi, Davies fc Alton 1999). 



Qcm(A) _ 2.005 (100 fim/Xy 
^) " 2390 h ~ (200 —< 



1/6- 



(12) 



Eqn.(|T^) has also been derived assuming that all grains 
have the same radius. The wavelength dependence of the 
ratio chan ges smoothly from A -1 to A~ 2 at 200 /im, as ob- 
served by Reach et al. (1995 ) using high S/N FIR spectra 
of the Galactic plane. The absolute value of the ratio has 
been derived using d ust column density maps calibrated to 
Galactic extinction ( Schlcgcl, Finkbcincr fe Davis 1998 ). 

Finally, FIR images are created integrating the emis- 
sion coefficient, 



B x (T d (r, z)),(13) 



j\(r,z) = n(r, z)a aba (X)Bx(T d (r, z)) 

= ~^k~ exp(-^/a<i- | z | //3 d ) 

along a given line of sight through the dust distribution, 
under the assumption that dust is optically thin to FIR 
radiation. Using the emissivity of Eqn. ( [l2| ) this is justified 
for any model with reasonable values of ry- 

As for the optical images, the far infrared images are 
produced in units of surface brightness. FIR images have 
the same extent and resolution of the optical images, i.e. 
a region of 12x12 stellar radial scalelengths around the 
centre of the galaxy is mapped in 201x201 pixels. 

5. A worked example: NGC 6946 

In the following we show an application of the radiative 
transfer and dust emission model to the spiral galaxy 
NGC 6946. The final goal is to find the parameters of 
the dust distribution compatible with the available obser- 
vations in the optical and FIR range. 

NGC 6946 is a large nearby Sc galaxy (A few basic 
properties are presented in Table ||) . The optical appear- 
ance of the galaxy is characterised by six prominent spi- 
ral arms of which the three arms originating from the 
NE quadrant are brighter and more developed than those 
in the SW ( [Tacconi fc Young 1990| ). [Trewhella (1998a|,bD 
finds that the NE interarm region suffers high extinction. 
The region appear fainter because of the dust effect, rather 
than being intrinsically less luminous. Indeed, polarised 
light is observed in the interarm regions as well as in the 
spiral arms (Fendt et al. 1998). 

NGC 6946 shows a centrally peaked molecular gas dis- 
tribution ( racconi fc Young 1989] ) dominant in the in- 



ner 10' over the flatter atomic gas component (Tacconi 
fc Young 1986| ). The galaxy is marginally resolved in FIR 



observations with KAO and IRAS (Engargiola 1991; Al- 



ton et al. 1998c), with the dust emission tending to follow 



the spiral arms and the bright central emission. The 117" 
resolution 200/im ISO image shows a morphology sim i- 
lar to the lOO^im IRAS observation (|Alton et al. 1998cj) . 



Higher resolution SCUBA images show a tight correlation 
between the dust emission at 850/xm and the dominant 



gas phase (Bianchi et al. 2000a). 

We now describe the observed SED for the stellar emis- 
sion, used in the normalization of the radiative transfer, as 



shown in Sect. 2.5, and the SED of dust emission, that will 



be compared to the model results. The procedure adopted 



in the modelling are presented in Sect. 5.3 



5.1. The stellar Spectral Energy Distribution 

We have derived the observed stellar SED from literature 
data. Values for the flux inside an aper ture of 5' (corre- 
sponding to the B-band half light radius; Engargiola 1991 ) 
are presented in Table [|. All the data have been corrected 
for Galactic extinction (Table |J) using the assumed ex- 
tinction law. 

Optical and Near Infrared data in t he bands U, B, g, V, 
r, I, J, H, K are from Engargiola (1991 ). The observed NIR 
emission is considered as purely stellar. Using the IRAS 
12 fim flux as a template of small grain emission and the 
Desert ct al. (1990| ) model, we derived a dust contribution 
to the K-band emission of only 0.3%. The stellar emission 
at 5 /^m has been extrapolated from the K-band, using the 
Rayleigh-Jeans spectral region in the synthetic galactic 
SEDs of |Fioc fc Rocca-Volmcrange (1997| ). 



The SED in the non-ionising UV is taken from Rifatto 



Longo & Capaccioli (1995b). The authors derive fluxes for 
a large sample of galaxies, by homogenising observation 
from several satellites (notably IUE), balloon and rocket- 
borne experiments, with different apertures and sensitiv- 
ities. An aperture correction depending on the morphol- 
ogycal type is applied to each flux, and total magnitudes 
are given for three photometric bands centred at 1650A 
(short-UV), 2500A (medium-UV), 3150A (long-VV). Us- 
in g an appropriate aperture correction and calibrating as 
in Rifatto, Longo fc Capaccioli (1995a ), we have derived 
fluxes for our selected aperture. The quoted errors come 
mainly from the aperture correction. The flux at the Ly- 
man limit (912A) has been extrapolated after observing 
that in a vf u versus log A plot the observed SED is flat for 
the short- and medium-\JV . This trend is assumed to be 
valid down to the ionization limit. The ionising UV is not 
included in the model (Sect. 0). 

5.2. The dust Spectral Energy Distribution 

The SED of dust emission is also shown in Table |[ Sur- 
face brightnesses over the half-light aperture shortward 
of 100/im have been derived from IRAS High Resolution 
(HiRes) images ( |Rice 1993|; Alton et al. 1998b|). De rived 
values have been colour corrected (Rice et al. 1988) and 



are consistent (within a 20% error; Alton et al. 1998c ) with 
the analogous data provided by Engargiola (1991 ), derived 
on previous enhanced resolution IRAS images. The value 
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RA (J2000) 


20 h 34 m 52:0 


Dec. 


+60° 9' 15" 


inclination 


34° 


P.A. 


64° 


distance 


5.5 Mpc 


D 25 


11.5' 


A B 


1.73 


Ldust 


1.5 - 3 -10 10 L 


M gas 


9 • 10 9 M 



Table 2. Basic properties of NGC 6946. 



de Vaucouleurs ct al. (1991), RC3 



RC3 



Garcia-Gomcz fc Athanassoula (1991 ) 



Garcia-Gomez fc Athanassoula (1991 ) 
Tully (1988D 



RC3 
RC3 

this work, 



Sect. 



Devereux & Young (199C), rescaled to 5.5Mpc 



stellar emission 



A 
/im 



C (5'©) 

Jy 



ohort UV 



0.091 
0.165 



0.05 
0.00 ± 0.02 



i edh 



JJM — 0.250 — 0.13 ± 0.09 



long-VV 
U 
B 

g 
V 
r 
I 
J 
H 
K 



0.315 
0.360 
0.435 
0.495 
0.554 
0.655 
0.850 
1.250 
1.650 
2.200 
5.000 



0.3 ± 0.2 
0.50 ± 0.05 
0.90 ± 0.06 

1.1 ± 0.1 
1.38 ± 0.07 

2.0 ± 0.2 

3.2 ± 0.3 

5.1 ± 0.4 

7.2 ± 0.6 
5.8 ± 0.5 

1.3 



dust emission 



12 


11 ± 2 


25 


14 ± 3 


60 


120 ± 20 


100 


240 ± 50 


160 


310 ± 60 


200 


280 ± 80 


450 


9-27 


850 


2-3 



Table 3. NGC 6946 mean flux inside a circular aperture of 
diameter 5' (corresponding to the B-band half light radius; 
Engargiola 1991), for several wavelengths. The origin of 
the data is described in Sect. 5.1 and 5.2, for the stellar 
and dust emission, respectively. 



at 160/iin has been derived by Engargiola (1991 ) from the 
air-borne KAO telescope. 

Data at 200/im are derived from P32 images taken 
from the ISOPHOT instrument aboard the ISO satellite 



(Alton et al. 1998c). The P32 mapping mode is still not 
scientifically validated and its photometric calibration is 



highly uncertain (Klaas & Richards 2000). Alton ct al.| 
(1998c|l compared both the integrated flux with the value 
derived by Engargiola (1991) on 200/im KAO images and 



the measured background with an extrapolation from the 
100/xm value. They conclude that the calibration may 



overestimate the flux by about 30%. A 30% error is shown 
in Table |. 

Finally, lower and upper limit of the flux at 450 and 



850 /im are derived from SCUBA images (Bianchi et al 



2000a). The lower limit correspond to the signal com- 
ing from regions 3<r brighter than the sky, inside the se- 
lected aperture. The upper limit is derived assuming a la 
emission for the regions without detected signal. Sub-mm 
fluxes are only given for completeness, the FIR models be- 
ing constrained mainly by the fluxes at 100, 160 and 200 
/im. Emission at 12, 25 and 60 /im, dominated by small 
grains (Desert ct al. 1990), is used in Sect. ^| to validate 
the MIR correction. 



5.3. Modelling procedure 

After choosing the star-dust geometry, a radiative trans- 
fer simulation is run for each of 17 bands of the model 
(Sect. |2.5| ) and images for the inclination of NGC 6946 
are produced (Table ||). An intrinsic scalelength a* — 2.5 
kpc has been used for the stellar emission, derived from 



a K-band image (Trcwhclla 1998b), where the effects of 
extinction are smaller. For the a*/ ((3+) ratio of Sect. |2.l| , 
this leads to f3± =170 pc. The simulations do not depend 
on the absolute values of the scalelengths, as long as they 
are shown as functions of scaled galactocentric and verti- 
cal distances. However, absolute values are needed if we 
want to derive correct values for the emitted and absorbed 
energies in each band. 

The simulated images are then integrated inside the 
half light radius aperture (derived from the B-band simu- 
lation) . Mean fluxes over each band of Table |l| are obtained 
by integrating a continuous SED derived from the data in 
Table 0. These fluxes are then used to scale the results 
of the aperture photometry on the simulation. Because 
of this normalisation, the stellar SED derived from the 
simulated images is the same for each model. The intrin- 
sic unextinguished energy emitted by stars is then derived 
from the Monte Carlo code. The half-light radius aperture 
was chosen in the earlier stages of our analysis because of 



the availability of one more UV data point in Engargiola 
|(199ip . However, if the distribution we have adopted is a 
correct description of the galaxy luminosity density, the 
normalization is independent of the aperture dimension, 
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being equivalent to assigning a value for po in Eqn. ([jj). 
After the normalization, FIR images for the selected in- 
clination are produced and the fluxes derived, as for the 
images of stellar emission. The simulated SED in the FIR 
is then compared to the observed one. 

For a sample of seven spirals, Alton et al 



(1998c) 



measured B-band and FIR scalelengths in a galactocen- 
tric distance range from 1.5' to 3.5' after smoothing the 
optical and the IRAS images to the poorest resolution of 
the 200/xm ISO map (FWHM=117"). For NGC 6946, they 
found that the B band scalelength is slightly smaller than 
the 200/Ltm scalelength, by a factor 0.9, while it is larger 
than the 100/im scalelength, by a factor 1.8. The abso- 
lute values of the scalelengths are presented in Table |]. 
The large 200/im scalelengths, compared to the optical 
and IRAS data, are a general property of the sample. The 
FIR simulation of this work are compared with the obser- 
vations of Alton et al. (1998c ). FIR images are smoothed 



to the ISO resolution and scalelengths are derived on the 
same distance range as the observations^. 

Keeping fixed the stellar distribution, several models 
are produced with different parameters for the dust dis- 
tribution. The final goal is to obtain a single model able 
to describe both the observed FIR SED and the spatial 
distribution of the emission. 

6. Results 

The first models presented here have the geometrical pa- 
rameters of the standard model, i.e. a star-dust geom etry 



2.2 



with ad = a* and /3d = 0.5/3*. As already said in Sect 
the choice is mainly motivated by the presence of extinc- 
tion lanes along the major axis of edge-on galaxies. The 
H2 column density in NGC 6946 has an exponential ra- 
dial profile, with a scalelength of 90" ( racconi fc Young 



1986). Incidentally, this is very close to a*. Thus, the dust 



disk discussed here can be thought of as a dust component 
associated with the molecular phase. Four values for the 
V-band face-on optical depth have been chosen, ry =0.5, 
1, 5 and 10. 

The SED of the four models in presented in Fig. 0. 
The thick solid line represents the stellar emission of the 
galaxy as derived from the simulations, normalised to the 
observed data points as described in the previous section. 
Each different line shows the intrinsic unextinguished SED 



in the short-wavelength part of the spectrum and the FIR 
emission in the long-wavelength side. For the dust-free 
emission, the SED is derived from the total intrinsic en- 
ergy, assuming isotropic emission and measuring the half 
light radius in a transparent model. All the models pre- 
sented in this paper show a spike at A «2000A in the un- 
extinguished SED. This is due to the extinction feature at 



2 As described in Sect. 2.4, a region of 12a* is covered by 
201 pixels, thus giving a pixel size of 5.7" (150 pc), for the 
adopted a* =95". Therefore the ISO beam can be modelled by 
a gaussian of FWHMf%20 pixels. 



1.5 



- 1 




o d =l.(.) a t 
(3 4 =0.5 (S. 



/ \ 



0.1 1 10 100 1000 

Fig. 2. Flux inside the B-band half-light radius for models 
with ay = a*, /3d = 0.5/3* and optical depths ry =0.5, L 
5 a nd 1 0. The data points are those described in Sect. 5.1 
and 5.2, for the stellar and dust emission, respectively. The 



solid line represents the SED of stellar emission measured 
on the simulated images. For each model, the lines in the 
UV-optical-NIR range represent the fluxes measured for 
a transparent model with the NGC 6946 inclination. Be- 
cause of the quantities plotted, vf v vs log(A), the area 
under the curve is proportional to the emitted energy. 



2175A present in the assumed extinction law (Sect. |2.3| ), 
while the observed SED is flat. A possible dip in the stellar 
SED caused by this extinction feature could be masked by 
the broad band of the observation or the large errors in 
the UV photometry. The extinction feature could also be 
absent, in which case our assumption of Galactic dust for 
NGC 6946 is not correct. 

In all the models, the wavebands from EUV to U con- 
tributes tho 30% of the total absorbed energy, while the 
wavebands from B to LMN to 60%. Therefore, most of the 
radiation is absorbed from the Optical- NIR light. A simi- 



lar result was obtained for the same galaxy by Trewhella 



(19981) 



The MIR corrections are quite similar for all the mod- 
els: approximately 32% of the total absorbed energy is 
estimated to go into MIR emission, the remaining 68% 
being available for thermal equilibrium processes and FIR 
emission. Since small-grains and PAHs responsible for no- 
equilibrium processes have a higher absorption efficiency 
at shorter wavelength, the contribution of absorption from 
Optical-NIR wavebands is higher after the MIR correc- 
tion. So the radiation originally emitted at A > 4000A 
contributes ~ 70% of the FIR emission. The same result 
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Fig. 3. Dust temperature map on 
a meridian plane for models with 
ay = a*, /3d = 0.5/3* and opti- 
cal depths Ty =0.5, 1, 5 and 10, 
from top to bottom. Temperature 
contours are plotted every IK and 
highlighted at 20 K and 30K by a 
label and a thicker line. The scale 
along the z-axis has been expanded 
for clarity. 



applies to all the models presented in this paper, therefore 
MIR corrections are not discussed separately in each case. 
A comparison between the estimated and observed MIR 
emission is given in Sect. 

The temperature distributions for each model are 
shown in Fig. |L as a function of the galactocentric radius 
and height above the plane. Apart from the central region, 
the distributions are very similar. For a dust distribution 
narrower than the stellar one, the stellar radiation field is 
expected to increase with height above the plane in an op- 



tically thick model (Draine & Lee 1984; Rowan- Robinson 
1986|)7|because the stars closer to the plane are shielded. 



This is evident in the central regions (R < 1.5a*) of the 
models. When the optical depth increases, dust at higher 
temperature is found at higher positions above the plane. 
In the models with higher extinction, the effect can still be 
seen at larger galactocentric distances, the region of higher 



temperature approaching the galactic plane at large dis- 
tances. Vertical gradients are very shallow, because of the 
greater extent of the galaxy in the radial direction with 
respect to the vertical and because the stellar distribution 
is smooth. 

The FIR spectrum is shown in Fig. ||. Only models 
with central optical depth between tv=5 and 10 produce 
enough energy to match the observational data. This is 
a general property of all the models we are going to dis- 
cuss: a substantial extinction is necessary to produce the 
observed SED in the FIR. The total amount of energy 
absorbed and re-emitted (in both MIR and FIR) is be- 
tween 1.5 and 3 • 10 10 L Q , for the two high Ty models, 
respectively. This corresponds to a fraction 0.3-0.4 of the 
intrinsic energy produced by the stars. Therefore, ^1/3 
of the bolometric luminosity of NGC 6946 is absorbed by 
dust. 
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ttB 


QlOO/jm 






observed (Alton et al. 1998c) 


1.52 


0.86 


1.66 




Standard model rv=0.5 


1.15 


0.65 


0.80 




Tv=1.0 


1.19 


0.66 


0.81 




t v =5.0 


1.63 


0.72 


0.86 




TV=10. 


2.07 


0.77 


0.91 




«d = 1.5a* model Ty—0.5 


1.16 


0.76 


1.01 


IS! 


7V=L0 


1.21 


0.76 


1.02 




r v =5.0 


1.60 


0.83 


1.08 




rv=10. 


1.76 


0.87 


1.12 




r v =0.5 HI +r v =5.0 H 2 distr. 


1.59 


0.81 


1.06 





Table 4. B-band, 100 fim and 200 fim scalelengths for 
several models, in units of a* (95"). Scalelengths have 
been measured on the models as described in Sect. 5.3. 



The unsmoothcd B-band scalelength, measured on a B 
band image ( [Trewhella 1998b| ), is 1.3a*. 



The optical (B) and FIR (100/mi and 200 fim) scale- 
lengths measured on standard models are shown in Ta- 
ble ||. The B-band scalelength increases with Ty and it is 
close to that observed if Ty ~ 5. The FIR scalelengths 
increase with the wavelength. A slight increase with Ty 
is also observed, because of the smaller temperature in 
the centre. However, FIR scalelengths are never as large 
as observed, especially at 200/im. While the 100/im scale- 
length for the ry = 5 model is underestimated by ~15%, 
the 200/im scalelength is about one half of that measured 
on ISO images. Therefore, the ratio between the B-band 
and 200/im scalelen gths is smaller than that measured by 
Alton et al. (1998c ). It is interesting to note that, because 
of the more rapid increase of the B-band scalelength with 
ry, the lowest values of the ratio are obtained for optically 
thin cases. Nevertheless, optically thin cases are unable to 
explain the large 200/im scalelength and the energy out- 
put. We note here that, by adopting a single a* measured 
in the K-band, we have assumed that the larger optical 
scalelength with respect to that measured in the NIR is 
entirely due to dust extinction. If indeed an intrinsic dif- 
ference exist, a smaller optical depth may be able to re- 
produce the observed B band scalelength. We will discuss 
the effect of our assumption in Sect. 0. 

Alton et al. (1998c) suggest that the large observed 
scalelengths at 200/im is due to a dust distribution more 
extended than the stellar disk. To test this hypothesis, we 
run models with ay = 1.5a* ( Davics et al. 1997|; Xilouris 
et al. 1999| ) , keeping the other parameters as for the stan- 



dard model. The ry = 5 case is shown in Fig. |I| As for 
standard models, only optically thick cases are able to 
match the observed SED. For the same optical depth the 
extended model has a higher extinction (e.g. 44% of the 
energy is absorbed in the V-band against the 34% of the 
standard model). 




= 5.0 standard 



T.-0.5 HI 
T v =5.0 H P distr. 



10 

X [/im] 



Fig. 4. Same as Fig. 0, but for the models with extended 
dust distributions. The SED for a standard model with 
the same optical depth is also included, for comparison. 



The temperature distribution for the extended model 
is shown in the central panel of Fig. ||. For ease of compar- 
ison, the temperature distribution of the standard model 
is shown again (in the right panel) , with the same scale as 
for the new model. Within a radius of 6 a* (the extent of 
the stellar disk) the temperature pattern of the extended 
model is quite similar, apart from a small difference due 
to the normalisation. This is reflected in the peak of the 
FIR SED, that is essentially the same in both the models. 
Outside of 6a*, the truncation of the stellar distribution, 
dust is colder and it does not modify the shape of the 
SED. The steeper temperature gradient at 6a* is an arti- 
fact due to the truncation of the stellar disk. A truncation 
is indeed suggested by counts of faint stellar sources in 
the Galaxy (Sect. 2.1). A few tests have been conducted 



with stellar distributions truncated at the same distance 
as the dust disks, to avoid having dust in regions without 
local stellar emission. The steeper gradient disappears and 
a larger distance is needed to reach the same temperature. 
However, changes are small, the general trend in the tem- 
perature distribution and in the FIR emission distribution 
are essentially the same. It is interesting to note that for 
R > 6a*, dust is colder on the plane than above, because 
starlight is seen through higher optical depths along the 
plane. 

The extended dust distribution causes an increase of 
the FIR scalelengths with respect to the standard model 
(Table ^). Scalelengths at 100/im are quite close to those 
observed (within a 3% for the ry = 5 case). However, 
the 200/im scalelengths are still underestimated, by at 
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least 30%. B band scalelengths are similar to those for 
the standard model (the effect of an extended dust dis- 
tribution being appreciable only for high optical depths). 
As in the standard model, despite the increase of the FIR 
scalelength with Ty, B/200/im scalelength ratios are close 
to those of Alton et al. (1998c[) only in the optically thin 



dard and extended dust distribution. FIR scalelengths are 
larger in extended models with respect to standard mod- 
els. The 100/im scalelength is very similar to that derived 
on IRAS images, but the 200/im scalelength is always 
smaller than that measured by ISO. It is not possible to 



case. 



The analysis of Galactic FIR emission of Davies et al 
(1997|)]suggests that the vertical scalelength of dust as well 
is larger than the stellar one, by a factor of two. Models 
that included a thicker dust disk, however, do not produce 
better results than those presented here ( Bianchi 1999| ). 

Larger FIR scalelengths can be produced by extend- 
ing further the dust distribution. Davies et al. (1999b) 
note that the IRAS emission tends to follow the molecular 
gas, while the 200/j,m profile is closer to the much broader 
atomic gas component. Could dust associated with the HI 
distribution be responsible of the emission at 200/im? A 
model was produced with two dust disks. As said earlier, 
a disk with ad = a* can describe dust associated with the 
molecular gas. As for the dust associated to the atomic 
gas, we chose a disk with a flat radial distribution up to 
2.5a*, then falling off exponentially with a scalelengt h of 3 
a*. This mimics the observed column density of HI (Tac- 
Young 1986| ). 



produce B/200/im scalelength ratios as small as in Alton 
et al. (1998cD , for any of the models presented here. 



com 



Using t he relation between extinction and gas column 
density of Bohlin et al. (1978) and the column density 



measured by Tacconi fc Young (1986 ), we derived an op- 
tical depth Ty ~ 0.5 for the broader dust disk. For the 
dust disk associated with the molecular phase we adopted 
Ty — 5, as in the previous standard disks. Using the 
observed H 2 column density, a value larger at least by 
a factor of two would have been derived. However, too 
much energy would have been absorbed from dust (see, 
for instance, the SED of the Ty = 10 standard model in 
Fig. ||). The discrepancy between the value derived from 
the gas column density and that needed to produce the 
right amount of absorption may be due to the smaller ef- 
fective absorption in a clumpy structure of dust. Clumping 
is discussed in Sect. [t]. 

The SED of this new model is shown in Fig. ||. The 
temperature distribution is presented in the bottom panel 
of Fig. ||. Despite the colder temperature of the dust asso- 
ciated to HI, the behaviour of the FIR radiation in the re- 
gion where scalelengths are measured is dominated by the 
dust disk associated with The derived scalelengths 
are similar to those of an extended model (ad = 1.5a*) 
with a single disk. 

In conclusion, the optically thick regime is required to 
match the observed FIR SED, both in models with stan- 



3 By extending the dust disk further, the FIR scalelength can 
became as large as that observed. This happens, for example, 
for a single Hl-like dust distribution. However, optically thick 
models (tv = 2 for the Hl-like distribution) are still needed. 
Such models would have a dust mass larger than what implied 
assuming the canonical gas-to-dust mass ratio (Sect. (?|). 



7. Discussion 

As shown in the previous sections, models with a central 
face-on optical depth Ty ~ 5 are necessary to explain the 
SED observed in the FIR for NGC 6946. A high optical 
depth through the central region s of the galaxy has also 
been f o und by Engargiola (1991 ) and Devereux fc Young 
(1993| ). [Evans (1992| ) and |Trewhella (1998bj ) apply the en 



ergy balance method to the stellar and dust emission of 
NGC 6946, using a TRIPLEX model, i.e. an analytical 
approximation for the radiative transfer, neglecting scat- 
tering, in a standard mode (Disney et al. 1989). They 



both derived high optical depths for the disk, using the 
data inside the half light radius. Evans (1992) measured 
ry = 6 - 7, while |Trewhella (1998b| ) Ty =~T± 1. A high 
optical depth is also suggested by high-resolution sub-mm 
SCUBA images: the diffuse inter-arm emission in the NE 
spiral arms at a distance of 2' (« a*) is compatible with 
Ty = 2.2 (|Bianchi et al. 2000a|). |Xu fc Buat (1995|) car- 



ried out an energy balance on a sample of 134 nearby 
spirals with available UV, B and IRAS fluxes. They de- 
rived a mean optical depth tb = 0.60. A direct compar- 
ison between their result and the optical depth needed 
by our model to explain the FIR output may lead to the 
wrong conclusion that we have overestimated the amount 
of absorbed energy. However, the mean ratio of bolomet- 
ric luminosity absorbed by dust for their sample (~ 1/3) 
is similar to that predicted by our models for NGC 6946. 
The difference depend on their adopted geometry, a plane- 
parallel homogeneous model for dust and stars. If the dust 
is associated with the dominant gas component, like obser- 
vations suggest, our choice of an exponential distribution 
seems more appropriate. For the same optical depth, their 
model results much more effective in extinguishing radia- 
tion. For example, in the UV, Ku & Buat adopt dust and 



stellar distributions of the same thickness (a slab model 



Disney et al. 1989) and isotropic scattering with albedo 



w=0.18. By opportunely setting scalelengths and trunca- 
tions, the BFG model can produce results for a slab. At 
their UV reference wavelength A = 2030A, 60% of the 
radiation is absorbed by dust, for their adopted UV ex- 
tinction law. Because of their low albedo, scattering does 
not play a relevant role in the radiative transfer. In our 
corresponding UV4 band, the Ty = 5 standard model ab- 
sorbs only 50% of the radiation. In the optical, they adopt 
a sandwich model, with the thickness of the stellar distri- 
bution twice that for dust. At A = 4400A, their model 
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Fig. 5. Same as Fig. ||, but for the extended models with ad = 1.5a* and ry =5 (central panel) and with two dust 
disk associated to each gas component (lower panel; see text for details). The temperature distribution for a standard 
model with ry =5 is shown in the top panel, for comparison. All the models have the same scale along the radial axis. 



absorbs 20% of the radiation, while the ry = 5 standard 
model 40%. 

In our simulations most of the radiation is absorbed 
in the Optical-NIR spectral range (60% for lig ht at A > 
4000 A). This was alredy note d for NGC 6946 by |Trcwhclla 
(1998b). [Xu fc Buat (1995| ) reach opposite conclusions, 
with the non-ionising UV (912A< A <3650A) contribut- 
ing to 60±9% of the absorbed radiation. This is in part 
due to the model differences we have already discussed. 
Xu & Buat adopted a different thickness for the UV stel- 



lar distribution to simulate the thinner distribution of 
younger stars. We tried a similar approach by adopting the 



same scaleheight for dust and star at smaller A (Bianchi 
199' ~ 



^|but the fraction of radiation absorbed in the UV 
increased only by a 5%. Furthermore, the SED observed 
in NGC 6946, with a peak in the NIR, may be different 



from the mean characteristic of the Xu & Buat sample 



The UV selected sample may be biased toward bright UV 



galaxies, although this hypothesis is dismissed in Buat & 
Xu (19f|). 

Another major difference between our model and that 



of Xu & Buat (1995) is in the ionising UV longward of 
912A. In their work the absorption of Lyman continuum 
photons contributes as much as 20±1% to the total FIR 
emission in a sample of 23 late type galaxies. Our omission 
of this spectral range from the radiative transfer, due to 



the lack of empyrical data, may result in an overestimate 
of the optical depth. To test our approximation, we de- 
r ived the ionising UV from Ha observations of NGC 6946 
( Kcnnicutt fc Kent 1983 ), assuming stan dard ionisation 
conditions as in HII regions ( [Lcqucux 198fj| ). After correct- 
ing for the [Nil] contamination (25% for spirals), Galactic 
extinction (An a ps 1) and internal extinction (30% for 
the R-band, where the Ha line is located, for a standard 
Ty — 5 model) , we derived an intrinsic unextinguished Ha 
flux, 



/(Ha) = 9.1 10 _11 erg cm- 2 s" 



(14) 



Following Xu & Buat (1995), the Lyman continuum flux 
can be derived as 



/(Lyc) = 33.9/(Ha) 



3.1 lO^erg enr 2 s" 1 , 



(15) 



of which 75% is assumed to be absorbed by gas and con- 
verted into emission lines at larger wavelengths (see also 



Mezger 197g; DeGioia-Eastwood 1992). Allowing the re- 



maining 25% to be entirely absorbed by dust and adopting 
a lower limit MIR correction of ss 70% (Sect. |]), the total 
FIR luminosity originating from the absorption of ionising 
photons is 



L FIR (Lyc) = 2.2 10 8 L e 



0; 



(16) 
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for the assumed distance of 5.5 Mpc. This correspond to 
only « 2% of the total FIR luminosity emitted by the stan- 
dard Ty—h model. As a comparison, the contribution to 
the FIR from the EUV band for the same model is 3.6%. 
The ionising UV contribution is similar for all the models 
that roughly provide the same amount of FIR energy as 
that observed. Therefore, it is justified to omit the ionising 



ISRF is larger than the local. As already said, emission 
from small grains can contribute to part of the observed 
FIR flux. However, for a wide range of heating conditions, 
the contribution close to the peak of FIR emission is min- 



imal Desert et al. (1990). Since modelled and observed 



radiation in the model. The approach of Xu & Buat (1995) 



is different from the one presented here. Their UV con- 
tribution includes direct absorption of Lyc photons and 
indirect (via emission lines), while in the present model 
the absorption of emission line photons is taken care of 
in the spectral band where the emission occurs (e.g. in 
the R-band for the Ha line), the total contribution of re- 
combination is summed up to the stellar SED for each 
band. Nevertheless, their ratio between Lyc emission and 
total absorbed energy is similar to the one derived here. 

A wrong MIR correction could also be responsible for 
an underestimation of the FIR emission in our model. The 
fraction of absorbed energy that goes into MIR emission 
depends essentially on the absorption of light from the 
short wavelength spectrum, since the absorption efficiency 
of small grains responsible for non-equilibrium processes 



fluxes are compared in this spectral range, the assump- 
tion that all FIR emission occurs at thermal equilibrium 
does not affect sensibly our results. 

The position of the peak emission in the modelled FIR 
SED clearly shows that the dust temperature in the sim- 
ulation s is not severely d ifferent from the real one. Re- 
cently, Stickel et al. (200C ) extracted high signal-to-noise 
sources from the ISOPHOT Serendipity Survey. Fluxes at 
170 fim were measured for 115 sources having a galaxy as- 
sociation, most of which are spirals, and integrated with 
100/im IRAS data. The distribution of the ratio of 170/im 
and 100/1 fluxes is quite well confined, with half of the 
galaxies having a ratio between 1 and 1.5. For the emissiv- 
ity adopted here, this translates in a colour temperature in 
the range 23-28 K. Similar colour temperatures (24-26K) 
can be derived from the ratios of the total fluxes at 100 
and 200/im, in any of our models. Using emissivities de- 



is high er in the UV (Sect. | D|). For the model presented 
here, with a duat acalchcight amallcr than the atcllar one, 
the amount of energy absorbed from UV bands does not 
increase very much with the optical depth (The saturation 
effect, Bianchi et al. 1996). The MIR corrections for these 
models are therefore quite constant, ~32% of the total 
absorbed energy being re-emitted in the MIR. Even for 
models with higher efficiency in extinguishing radiation, 
like the optically thick model with a thicker dust distribu- 
tion, the fraction of the total energy emitted in th e MIR 
is not very different from this value ( Bianchi 1999 ). 

For the local ISRF the contribution of small grain cmis- 



rived from Galactic emission and exinction, Bianchi et al 
(1999|) found that the value of the dust mass for a small 



sample of galaxies does not depend dramatically on the 
assumed spectral behaviour of the emissivity. This can be 
easily interpred if dust temperatures in external galaxies 
are similar to those in the Milky Way. Because of these 
similarities, and of the small dependence of the temper- 
ature on the assumed model, it is not too audacious to 
compare values measured on the present models with the 
high quality observations of the Galaxy. 

Dust temper atures at high latitude are quite constant 
in the Galaxy (Reach et al. 1995; Lagache et al. 1998 



sion to the 60 /im IRAS band is ^62%, while at 100 /im Schlcgcl et al. 1998| ). For the emissivity adopted here, 



it is o nly 14% and at 200 /im 4% (Desert et al. 1990|). a tem perature of fa21K can be derived ( Bianchi et al 



Theref )re, the fraction of energy emitted in non equilib- 
rium heating can be roughly estimated by measuring the 
MIR emission shortward of 60 /im. After integrating a con- 
tinuous SED interpolated from the data points in Table ||, 
the MIR energy is derived to be 34% of the total infrared 
energy emitted by dust. The value derived from obser- 
vation is very close to the model one. This justifies the 
use of the Desert et al. (1990| ) dust model as described in 
Sect, g It is interesting to note that the infrared galactic 
spectrum used in the Desert et al. ( 1990| ) model is differ- 
ent from the one of NGC 6946. As an example, the ratio 
between fluxes at 60 /tm and 100 /im is 0.2, while it is 
0.5 from the NGC 6946 data. This does not necessarily 
mean that the dust model of Desert et al. (199fj| ) cannot 
be applied to NGC 6946. The different ratio could be due 
to different heating conditions in the local interstellar ra- 
diation field, with respect to the mean radiation field of 
NGC 6946. Larger r atios between 60 /jm and 100 /jm can 
be derived from the Desert et al. (199C ) model when the 



1999). Remarkably, nearly all the models present a sim- 



ilar temperature at a galactocentric distance of 3a* (the 
Sun posit ion, for i? Q =8.5kpc a nd a Galactic scalclcngth 
of 3kpc). [godroski et al. (1997 ) decompose the Galactic 
FIR emission observed by DIRBE into three components, 
associated with the atomic, molecular and the ionised gas 
phases, and derive the temperature for four annulii at dif- 
ferent galactocentric distance. Fig. ^ shows the temper- 
ature of the dust associated with the atomic gas (sup- 
posedly heated mainly by the diffuse ISRF) for each an- 
nulus. Data have been scaled to a Galactic scalclcngth 
a* = 3 kpc and corrected for the emissivity law used in 
this work. In Fig. |^ we also plot the temperature radial 
profile for four representative models: the Ty = 0.5 and 
Ty = 5 standard models; the Ty = 5 ad = 1.5a* model; 
and the model with a dust disk associated with each of the 
two gas phases. Clearly, all the simulations presented here 
have a temperature gradient compatible with the Galactic 
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v = 0.5 standard model 
, = 5.0 

*— 5.0 £X d — 1.5a. model 

■ =0.5 HI + t =5.0 H, dist. 



standard models. Alton et al. derived a ratio 0.9 between 



H [a.] 



Fig. R.I Temperature along the gaJactk plane as a. function 



of the galactocentric distance for four different models: the 
Ty = 0.5 and Ty = 5 standard models, a Ty = 5 «d = 
1.5a* model and the model with a dust disk associated 
with each of the two gas phases (Sect. |6|). Data points 
show the radial gradient of temperature in the Galaxy 
and are derived from Bodroski et al. (1997| ) as described 
in the text. The vertical error bar represent the error in the 
temperature determination, while the horizontal the radial 
range for which the temperature has been measured. The 
cross marks the temperature of 21K at the Sun distance 



from the galactic centre (Bianchi et al. 1999) 



Because of the lack of high resolution observations in 
the FIR, the derivation of temperature gradients in exter- 
nal galaxies is more difficult. Davies et al. (1999b| ) derived 
the temperatures from 100/im IRAS and 200/im ISO fluxes 
at two different position on NGC 6946, in the centre and 
on the disk at 3 arcmin from the centre. After smoothing 
for the ISO resolution, our models are compatible with 
those observations ( Bianchi 1999| ). 

For any model, the FIR scalelengths increase with the 
wavelength of emission, the ad/a* ratio and, although 
slightly with the optical depth. Therefore, larger FIR 
scalelengths can be found in optically thick extended mod- 
els. For the Ty = 5 extended model, that provides a good 
fit to the SED, the 100/mi scalelenght is very close to the 
value derived from IRAS images. The 200/xm scalelength 
is larger than scalelength at 100/im. However, it is 30% 
smaller than that observed by ISO. It has not been pos- 
sible to find a model able to reproduce the large 200/im 



scalelength measured by Alton et al. (1998c) in a sample 



the B and 200/im scalelengths. Because of the dependence 
of the B-band scalelength on Ty, similar ratios can be ob- 
tained only for optically thin extended models. Such mod- 
els fail to predict both the SED and the absolute values 
of the FIR scalelengths. A possible reason for the discrep- 
ancy between observations and models may reside in the 
transition effects of the ISOPHOT detectors. Memory ef- 
fects during the scanning of a bright source affects the P32 
mapping mode. Alton et al. checked for this problem by 
conducting the scalelength analysis on both in-scan and 
cross-scan directions and concluded that the effects were 
negligible. However, a proper analysis requires a knowl- 
edge of the underlying light distribution and a description 
of the transient effects, which are still poorly understood. 
The P32 mapping mode is still not scientifically validated 
( |Klaas fc Richards 2000| ). 

We have carried out a few test to check for the in- 
fluence of different stellar distributions, like including a 
small bulge appropriate for a late type galaxy or having 
a smaller vertical scalelength for younger stars (Bianchi 
1999). The result were not significantly different from 



of seven galaxies including NGC 6946. The B-band ratio 
is close to that observed for the Ty = 5 extended and 



those presented here. A basic assumption of all our models 
is that the intrinsic (not extinguished) radial stellar scale- 
length is the same at any wavelength. According to this 
view, the larger scalelengths observed in the optical are 
due to extinction, rather than to intrinsic color gradients. 
As outlined in Sect. 2.1, the problem is still very debated. 
A simple test can be conducted here to assess the influ- 
ence of scalelength variations with A. Let's assume that 
the measured B-band scalelength (1.3a*) is more repre- 
sentative of the intrinsic radial variation of the stellar 
distribution. By substituting a* with 1.3a*, the results 
of Table |Jshow that standard or extended models with 
Ty = 0.5 or 1 have an observed (after smoothing) B-ban d 
scalelength quite similar to the one of Alton et al. (1998c ), 
the difference between observed and intrinsic scalelength 
due to the smoothing process, rather than to extinction. 
Since the modelled B /200/im scalelength ratio remains the 
same, these optically thin cases may have a 200/im scale- 
length close to the observed. However, the fraction of en- 
ergy absorbed in these models will remain the same as 
well. Therefore, the optically thin models will still be un- 
able to produce the required energy output in the FIR. 
Furthermore, most of the energy is absorbed in the NIR, 
where the peak of stellar emission occurs. In this spectral 
range our K-band based a* is surely a better description 
for the galaxy radial scalelength. 

We also produced a model with two dust distributions, 
a ry = 5 standard disk associated with the molecular gas 
dominant in the central part of NGC 6946, and a Ty = 0.5 
disk associated with the broader distribution of atomic 
gas. Even with the presence of such an extended distribu- 
tion, the double disk model is dominated by the optically 
thick disk associated to the H2 necessary for the FIR emis- 
sion. For what concerns FIR scalelengths, the same results 
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as for the Ty = 5 extended model are obtained. It is in- 
teresting to note that for the parameters adopted here for 
the two dust disks and assuming the atomic + molecu- 
lar gas mass of Table |[ the gas-to-dust mass ratio is 180, 
close to the Galactic value of 160 (sodroski et al. 1994). 
Similar values can be retrieved from all the optically thick 
models. 

The high optical depth of the models contrasts with 
the recent determination of optical depth in edge-on spi- 



rals by pCilouris et al. (199S ). Using a sample of seven edge- 
on galaxies, they find a mean central face-on optical depth 
ry ~ 0.5. The higher opacity of NGC 6946 may be a result 
of the galaxy being very gas- rich ( Tacconi fc Young 1990] ) ; 



or it may be due to clumping of the ISM, affecting in a 
different way FIR and optical determinations of the opti- 
cal depth. While FIR observations detect all of the dust 
(at least when the temperature of the clump and inter- 
clump medium are similar), optical observations may be 
affected preferentially by the extinction of the smoother, 
lower density (and optical depth) inter-clump medium. 

To our knowledge, only two works in literature discuss 
the effects of clumping within the framework of a radiative 
transfer model for spiral galaxies. Kuchinski et al. (1998 ) 
distribute clumps in the dust disk assuming a constant 
filling factor and a value for the ratio between densities 
in clumps and in the nearby smooth medium. The other 
model is by Bianchi et al. (2000b), based on the BFG ra- 
diative transfer code, as the present work. Bianchi et al 



(2000b) modelled the clump distribution in a similar way 
to the distribution of molecular gas in the Galaxy. The 
clump properties were derived from those observed in Gi- 



ant Molecular Clouds. A comparison with the Kuchinski 
et al. ( 1998| ) results indicates a strong dependence of the 



observed brightness profiles on the detailed internal and 
spatial distribution properties of clumps. This makes the 
interpretation of the data very difficult. For the same dust 
mass, it is found that a clumpy dust medium has lower ex- 
tinction. However, when a fraction of the stellar radiation 
is allowed to be emitted from inside the clumps, as for Gi- 
ant Molecular Clouds hosting star-formation, extinction 
increases and can reach higher values than those for ho- 
mogeneous models. 

Predictions of the influence of clumping on the FIR 
models are not easy. If the dust component associated with 
the H 2 is clumped, it could be responsible for most of the 
FIR emission, when embedded stellar emission is consid- 
ered. A diffuse component associated with the HI would be 
responsible for a reduced Ty derived from the analysis of 
edge-on galaxies. For the models of Bianchi et al. (2000b| ), 



it is unlikely that a clumpy distribution with a large dust 
mass (corresponding to an optically thick smooth distri- 
bution) has an optically thin apparent ry as measured by 



Xilouris et al. (1999). However, the models of Bianchi et al 



(2000b) are based on the Galaxy, while clumping based on 
the dif ferent distributions of atomic and molecular gas in 
NGC f^946 may have a different behaviour. Clumping will 



also affect the spatial distribution of the FIR emission. 
Dust in clumps is shielded from the ISRF and heated to 
lower temperature. Depending on the distribution of cold 
dust clouds, a broader 200/xm emission could be produced. 

Two recent models include clumping of dust and em- 
bedded stellar emission to describe the radiative trans- 
fer and FIR emission of NGC 6946. However, the com- 
plexity of the models prevents an isolation of the effect 



of the dust distribution on the FIR heating. Silva et al 



(1998) fit the observed optical and FIR SED of NGC 6946 



with their galactic photometric evolution model, adopt- 
ing a simplified treatment for the radiative transfer in 
the diffuse medium and a separate model for dust in 
clumps. They find that nearly half of the dust emission 
comes from within molecular clouds. The adopted dis- 
tributions for smooth dust and stars are different from 
those in Sect. |[ the radial scalelengths are the same for 
both dust and stars, as well as the vertical scalelengths; 
the radial scalelength is twice the one adopted here; the 
disk is thicker by a factor of two. Since we use the same 
SED and their model is optically thin to B-band radiation 
(only 10% of radiation absorbed), most of the FIR emis- 
sion must come from absorption of UV photons in the 
molecular cloud distribution. A distribution of molecular 
clouds compatible with the galactic gravitational poten- 
tial, together with a smooth phase modelled on the atomic 
gas, is used by |5auty et al. (1998 ) to study the FIR emis- 
sion in NGC 6946. The radiative transfer is carried out 
using the Monte Carlo technique only for the UV radia- 
tion emitted in star-forming regions within the molecular 
clouds. The ISRF at A > 2000A is derived from a R-band 
map of the galaxy, scaled on the Galactic local ISRF. As 
in our case, the observed FIR emission can be explained 
only with substantial extinction. However, they find that 
emission at A < 2000A contributes 72% of the total FIR, 
while in the present work energy absorbed from optical 
radiation is dominant. 



Therefore, accordin g to the models of Silva et al. (1998 ) 
and Sauty et al. (1998] ) , the dominant contribution to the 
dust heating comes from young stars. If this is true, the 
recent star-formation rate of a galaxy could be measured 
using FIR fluxes. These are more readily available than 
other tracers of star- formation (Devereux & Young 1991). 
In a series of papers, Devereux & Young (1993, and refer- 



ences therein) compared FIR and Ha fluxes and suggested 
that the FIR luminosity is dominated by warm dust ab- 
sorbing radiation from OB stars. However, this result is de- 



bated. Walterbos & Greenawalt (1996) modelled the FIR 
emission in spiral galaxies, deriving the diffuse ISRF from 
optical profiles and estimating the dust column density 
from the atomic gas. For a sample of 20 galaxies, they 
found that dust heated by a diffuse ISRF can account 
for, on average, half of the observed IRAS fluxes. Using 
large nearby objects it is possible to decomp ose the FIR 
emission into the different sources of heating. Ku fc Helou 



(1996a ) measure the ratio between the IRAS 60fim and 
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Ha fluxes from bright FIR-resolved sources in M31. Using 
the total Ha luminosity and the Desert et al. (1990) dust 
model, they extrapolate the fraction of the total luminos- 
ity emitted by dust that is associated with HII regions 
and star-formation. A value of 30±14% only is found. A 
decomposition of the Galactic flux observed by COBE at 
140/im and 240/im shows that most of the dust emission 
(~70%) arises from dust associated with the atomic gas 
( [Bodroski et al. 1994 ), with temperature gradients com- 
patible with a diffuse ISRF heating. Only 20% of the FIR 
is emitted by dust associated with the molecular compo- 
nent and 10% is due to hot dust associated with the HII 
phase. 

The smooth models presented here attempt to explain 
the FIR emission without invoking dust heated in star- 
forming regions. Under the assumption of a diffuse stellar 
emission and dust, we find that radiation at A < 4000 A 
contributes « 30% to the total FIR, this including the 
ionising UV. This percentage will surely increase when a 
clumpy dust structure is adopted. The high optical depths 
necessary to produce the observed energy may be an over- 
estimation because we neglect clumpy hot dust. However, 
without a proper model for the radiative transfer and dust 
emission in a clumpy medium, it is not easy to under- 
stand if the results shown here are severely biased by the 
assumption of smooth distributions. 

Finally, we tried to model the FIR emission with an 
extra dust component, a spherical halo. The presence of 
dust in the halo of normal, non-active, spirals is predicted 
as a result of the imbalance between the radiation pressure 



and the galactic gravitational force (Davies et al. 1999a). 
A halo of cold gas clouds, able to explain part of the dark 
matter in spiral galaxies, may be stabilized by the pres- 
ence of dust grains ( pcrhard fc Silk 1996p . Unfortunately, 
it is difficult to obtain information about a putative dusty 
halo. Because such a distribution would act as a screen 



for the galactic disk, the radiative transfer fits of Xilouris 
et al. 1999) are unable to detect it. Analysing the dif- 



ference in colours of background objects between fields at 
different distan ces from the centre of two nearby galaxies, 



Zaritsky (1994 ) find a B-I colour excess in the inner fields. 



He derived a dust halo scalelength of 31±8 kpc, although 
the measured colour difference is only 2a of the statistical 
noise. 

Because of the high uncertainty of the halo parameters, 
we simply used a spherical homogeneous dust halo, to- 
gether with a standard disk. The halo has the same radial 
dimension as the dust and stellar disks, i.e. 6a*. Among 
the models we tried, the case in which Ty = 1 for both halo 
and disk provides a reasonable fit to the FIR emitted en- 



ergy and scalelengths (Bianchi 1999). However, when the 
galaxy is seen edge-on, such a dust distribution would be 
easily observable, within the resolution and sensitivity of 
available FIR observations. Alton et al. (1998b| ) studied 
the FIR emission in 24 edge-on galaxies, including star- 
burst and quiescent objects, using HiRes IRAS images. 



None of the object was found to be resolved along the 
minor axis. 



8. Summary and conclusions 

We have described in this paper a model for dust emis- 
sion in spiral galaxies, based on the Monte Carlo radiative 
transfer code of Bianchi, Ferrara & Giovanardi (1996). For 
each relative star-dust geometry and dust optical depth, 
the radiative transfer code is carried out for 17 different 
photometric bands, covering the spectral range of stellar 
emission. In the application to NGC 6946 presented here, 
we have chosen a model stellar SED and dust distribution 
that produces the same fluxes as those observed in the UV- 
Optical-NIR images. The code also produces a map of the 
total energy absorbed by dust. For each position inside the 
galaxy, dust is heated by an ISRF consistent with the ra- 
diative transfer itself, without any other assumption. The 
dust temperature is computed f rom the absorbed ene rgy, 
using the emissivity derived by Bianchi et al. (1999 ) for 
Galactic dust and correcting for the contribution of non- 
thermal equilibrium processes to the emission. Hence, FIR 
maps can be easily obtained for any wavelength, integrat- 
ing along a specific line of sight. 

The model optical and FIR scalelengths and the SED 
have been compared to NGC 6946 observations. Several 
models have been explored. It is found that optically thick 
dust disks (ry ~5) are needed to match the observed FIR 
output. Approximately one third of the total stellar radi- 
ation is absorbed by dust in this case. The temperature 
distributions are quite similar, for any of the dust disk 
models. Temperature values in the models are compara- 
ble with those observed in our Galaxy and other spirals. 

We have compar ed the modelled FIR scalelengths with 
the observations of |Alton et al. (1998c ). Extended dust 



disk model with ad = 1.5a* ( Davies et al. 1997; Xilouris 
|ct al. 19*99 ) have larger FIR scalelengths than standard 
models with ad = a*. For optically thick cases, the 100/im 
scalelength is close to the value derived on IRAS images. 
Alton et al. (1998c ) found that the B-band scalelength of 
NGC 6946 is smaller than the 200/xm one by a factor 0.9 
(0.8 for a sample of seven galaxies). We have not been 
able to reproduce the large FIR scalelengths measured on 
200/im ISO data. In the required optically thick regime, 
the scalelength ratio B/200/im is always larger than ob- 
served. Smaller ratios can be obtained only in optically 
thin cases, but the absolute values for the scalelengths are 
smaller. The results are not improved if two dust disks 
modelled on the gas distribution are used. The behaviour 
of the model is dominated by the standard optically thick 
disk associated with the molecular component, rather than 
the very extended dust distribution associated with HI. A 
spherical dust halo could produce results closed to those 
observed, but would also be easily detected in currently 
available FIR observations, which is not. 
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The high optical depth found for NGC 6946 con- 
trasts with recent determinations on edge-on spiral galax- 
ies Xilouris et al. (1999). This may be a result of our 
assumption of a smooth distribution for stars and dust. 
The inclusion of clumping in a proper model of radiative 
transfer and FIR emission is therefore desirable. However, 
the heavy dependence of clumping on the assumed model 
makes the modelling more complicate. Future high resolu- 
tion and sensitivity instrumentation will therefore be es- 
sential to dehne the dust distribution and limit the num- 
ber of parameters in the model. Clumping may also be 
responsible for the discrepancy between the observed and 
modelled scalelengths at 200 /im, if a diffuse component 
of cold dust clumps shielded from the ISRF is present. On 
the other hand, it is necessary to remind here that the 
results of |Alton et al. (1998c ) at 200/im are based on data 



which is not yet scientifically validated. Again, future FIR 
instrumentation or a set of validated data will help to asses 
if the large scalelengths are an artefact of the ISOPHOT 
detector transients or if they are the genuine results of a 
temperature distribution different from the one we have 
derived here. 

The main purpose of the paper was the presentation of 
the model itself. A test has been conducted on a well stud- 



ied galaxy. However, NGC 6946 is very gas rich (Tacconi & 
Young 1990| ) and its characteristics may be different from 



those of a 'mean' galaxy. Optical and FIR data are avail- 
able for several galaxies. A future paper will be devoted 
to their analysis and more general conclusions about the 
dust distribution and extinction will be drawn. 
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